Metallization of Nanofilms in Strong Adiabatic Electric Fields 
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We introduce an effect of metallization of dielectric nanofilms by strong, adiabatically varying 
electric fields. The metallization causes optical properties of a dielectric film to become similar to 
those of a plasmonic metal (strong absorption and negative permittivity at low optical frequencies). 
The is a quantum effect, which is exponentially size-dependent, occurring at fields on the order of 
0.1 V/A and pulse durations ranging from ~ 1 fs to ~ 10 ns for film thickness 3 — 10 nm. 



PACS numbers: 73.20. Mf 77.22.Jp 42.65.Re, 72.20.Ht 

Effects of strong electric fields on electron states in 
crystals have attracted a great deal of attention over 
many decades going back to Zener who predicted break- 
down due to interband tunneling 1 . In insulators this re- 
quires electric fields on the order of atomic fields £ ~ 
1 — 10 V/A. Interest to strong-field condensed matter 
physics has recently greatly increased due to the avail- 
ability of such strong electric fields in laser pulses of in- 
tensities I ~ 10 13 — 10 15 W/cm 2 . Ultrashort laser pulses 
with a few optical oscillations 2 * 3 - open up a possibility 
to study ultrastrong field phenomena in solids during 
periods of time too short for the lattice ions to move 
significantly. Recent ab initio calculations^ have repro- 
duced the Zener breakdown in insulators induced by a 
laser pulse of intensity ~ 10 15 W/cm 2 . Other strong- 
field phenomena that can be observable in crystals at a 
comparable field strength are the appearance of local- 
ized electron states, Wannier-Stark ladder in the energy 
spectrum^, and Bloch oscillations X At orders of magni- 
tude lower intensities, low- frequency optical fields cause a 
reduction of the band gap in semiconductors and insula- 
tors (Franz-Keldysh effect, FKE)^. The quantum con- 
fined FKE takes place in semiconductor quantum wells 
and is determined not by the field but by the total po- 
tential droppi£ It requires typical fields £ ~ 10~ 3 V/A. 

In this Letter we introduce an effect of metallization 
in insulator nanofilms, which is predicted to occur in ap- 
plied electric fields £ ~ 0.1 V/A. It is based on adia- 
batic electron transfer in space across the nanofilm. The 
minimum duration of the field pulse required for the adi- 
abaticity exponentially depends on the crystal thickness 
varying from ~ 1 fs for a 3 nm film to ~ 10 ns for a 
10 nm film thickness. This metallization effect manifests 
itself by a dramatic change in the optical properties of 
the system, which start to remind those of metals. In 
particular, plasmonic phenomena emerge. 

To demonstrate the metallization effect, we need to 
solve the one-electron Schrodinger equation for a periodic 
potential plus a uniform electric field very accurately. We 
will employ the widely used Kronig-Penney model for 
electrons in a film confined in the x direction by an infi- 
nite potential well. The corresponding potential energy 



(neglecting the electron-electron interaction) is 



V(r) 




U{y) + U{z) + e£x \x\ < L/2 
\x\ >L/2 



(1) 



where U(x) = — ot^ n= _ oa 5(x — na), and a is the lattice 
constant. Crystal thickness L is determined by the num- 
ber of the lattice periods N in the x direction, L = Na. 
Though this model does not precisely correspond to any 
real system, it is exactly solvable and catches the qual- 
itative features of the strong-field phenomena from for- 
mation of the quantum bouncer (QB) states to band gap 
collapse and metallization transition. We use the trans- 
fer matrix to find an exact solution of the Schrodinger 
equation with potential ([T]) - see Sec. II of Supplemental 
Material 11 . The zero boundary conditions for the wave 
function at z = ±L/2 have been imposed and energies 
found semi-analytically using the shooting method. 

The energy bands of the infinite crystal in the zero 
field are shown in Fig. [T] (a) by color coded lines for the 
three lowest bands: localized (green), valence (red), and 
conduction (blue). The discrete electronic levels for the 
N = 30 film in zero field are represented by the black dots 
superimposed on these infinite-crystal dispersion curves. 
Our computations are made for a — 13.5 eV A and 
a = 2.7 A to result in the gap between the conduc- 
tion and valence bands to be E g = 9 eV, which is the 
same as for Si02 . We consider an insulator (semiconduc- 
tor) film where the localized and valence bands are com- 
pletely filled, and the conduction band is empty. This 
corresponds to a population of 16 electrons per unit cell 
(taking spin into account) and a reasonable electron den- 
sity of n = 8 ■ 10 23 cm~ 3 . 

The energy spectrum of the nanofilm as a function of 
the applied field £ normal to the film is shown in Fig. Q] 
(b) where the filled bands are coded by the green (local- 
ized band) and red (valence band), and the empty con- 
duction band is indicated by blue. With an increase of £ , 
the linear Stark effect takes place, and the conduction- 
valence band gap decreases, completely closing at a met- 
allization field £ m = 0.148 V/A indicated by the vertical 
dashed line. This behavior can be understood from an- 
alytical theory presented in Sec. V of the Supplemental 
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FIG. 1: Field effect on electron states: spectrum, localiza- 
tion, transmission, and mixing. The data are for an N — 30 
(9 nm) nanofilm thickness, (a) Energy bands of infinite crys- 
tal (lines) and nanofilm (dots) in zero electric field. Bands are 
color coded by the green (localized band), red (valence band), 
and blue (conduction band), (b) Energy bands of nanofilm as 
a function of the applied electric field (color coded as above), 
(c) Valence band levels in electric field £ = 0.06 V/A. The 
black lines and state notations are discussed in the text, (d) 
The same as in panel (c) but for £ = 0.1 V/A. (e) Transmis- 
sion coefficient T as a function of electron energy E for dif- 
ferent applied fields (color coded as shown), (f) Anticrossing 
between the lowest-energy level in the conduction band (blue) 
and the highest-energy level of the valence band (red). Energy 
AE and field A£ are given with respect to the crossing point. 
(g) The wave functions of the anticrossing states of panel (f ) 
with the corresponding color coding for AS = 0.14 V/cm. (h) 
The minimum splitting of the anticrossing levels as a function 
of the film thickness expressed as N. The black dots are ob- 
tained from numerical computation. The red line is calculated 
using Eq. Q and the dashed blue line from Eq. ([5]). 



Material^, where Eq. (38) expresses £ m as 



E„ 



e[L - CiA e (S m ) - ^A h (S m )} 



(2) 



Here an electric-field quantum confinement length is 
A e ,fc(£) = [h 2 /(2m* eh e£ )] 1/3 , where to* h are the effec- 
tive masses for electrons and holes, and £i and Ci are the 
first roots of Ai(— x) and Ai (— x), respectively. 

The band gap can be estimated as E g ~ 
7T 2 H 2 J (2ma 2 ) . From this, we can estimate A e ,/, ~ 

1 1/3 

to/(7t 2 to* h)a 2 L <C L. Neglecting A e< h in compari- 



son with L, one obtains a very good approximation for 
£ m , band edges E^^ (where b and t stand for the top and 
bottom), and band gap E g as 

£m = -£ , E b}t (£) = £ M T e£- , E g {£) =E g - e£L . 

This implies a linear Stark effect near the metallization 
point, in an excellent agreement with Fig. [1] (b). 

An applied normal electric field causes the appearance 
of states localized between the corresponding boundary 
at X = ±L/2 and the stopping points whose coordinates 
x s and energy E are related as E — E^t + e£x s [ci. Eqs. 
(29) and (33) of the Supplemental Material^] . These re- 
lations are represented by the slanted black lines in Figs. 
Q] (c) and (d). Carriers at the conduction band bottom 
behave as electrons, and those at the top of the valence 
band behave as holes. The corresponding mobility edges 
are denoted as E\ and E 2 . For moderate fields [panel 
(c)], the states with energies E > E 2 or E < E\ are 
Bloch-electron quantum bouncers (QBs)i2. The states 
with intermediate energies E 2 > E > E\ do not have 
stopping points and are delocalized, conducting (DC). 

For a stronger electric field, as shown in Fig. Q] (d), 
the mobility edges overlap, E\ > E 2 , and the DC states 
disappear. Instead, localized Wannier-Stark (WS) states 
appear with energies E\ > E > E 2 . These states are very 
close to those in infinite lattices^ Their contribution to 
the static conductivity vanishes sec the next paragraph. 

Assume that the barriers at the boundaries of the 
nanofilm are transparent enough to allow for tunneling 
through. Then it is physically meaningful to find the 
transmission coefficient T for carriers injected at a cer- 
tain energy, which is plotted in Fig. [1] (e) for the valence 
and conduction bands for different values of £ . As £ in- 
creases, the transmission band collapses along with the 
disappearance of the delocalized conducting (DC) states. 

Metallic behavior is actually characterized by two dif- 
ferent phenomena: dc conductivity and negative Ree 
(where e is permittivity) that contributes to plasmonic 
phenomena. They do not necessarily both take place 
simultaneously. We have shown above that the field- 
induced localization eliminates the dc conductivity. 

Now we will show that, to the opposite, the applied 
field that increases slowly (adiabatically) turns an insu- 
lator crystal into a metal optically, which we call the 
metallization effect. A distinct property of metals is the 
absence of a band gap around the Fermi energy. The 
two levels at the edges of the valence and conduction 
band at field £ m experience an anticrossing as Fig. [1] (f ) 
shows. The band gap is very small, AE m ~ 0.1 /ueV, 
which stems from a very little overlap between the wave 
functions of the two edge state QBs [Fig. [T] (g)]. 

The minimum splitting of these anticrossing levels is 
related to the matrix element of the Zener-type tunneling 
between the valence and conduction bands. This splitting 
can be calculated analytically using the quasiclassical ap- 
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FIG. 2: Electron states in metallizing fields for N — 30. The 
occupied states (below Fermi energy Ef) are indicated by red 
and the vacant ones by blue, (a) Electrons in electric field £ = 
0.144 V/A slightly less than the metallization threshold S m — 
0.148 V/A. (b) Electrons in electric field £ = 0.152 V/A 
slightly exceeding £ m . The black double arrows denote the 
dominant low-frequency optical transitions, (c) and (d) The 
same as for previous panels but for fields £ = 0.238 V/A and 
£ = 0.854 V/A, correspondingly. 



proximation of Refs. [p - fla and Eq. ([3]). We obtain 



A£L 



ex v/f^cxp -7r J [lE* I '(AheSm) (4) 



ex exp 



-7rVA-V7(4ft) 



(5) 



where the reduced mass is /x = TO*m^/(m* + to£). 

The dependence of the minimum band splitting A.E m 
on N = L/a from Eq. (QJ, where the critical field is given 
by Eq. ([2]), is shown in Fig. [T](h) with a solid red line. It 
is in an excellent agreement with numerically computed 
points obtained from our quantum-mechanical solution, 
which are displayed as bold dots. Dependence given by 
a simplified expression ([5]) is shown by the dash blue line 
in Fig. [T] (h) and is a good approximation. 

Now we turn to the metallization. Consider first 
field £ = 0.144 V/A, which is slightly less than £ m = 
0.148 V/A. The corresponding energy levels near the 
Fermi energy Ep sa 14.8 eV are shown in Fig. [2] (a). 
We assume temperature to be sufficiently low so only the 
states below Ep are occupied (shown by red) and those 
with E > Ep are vacant (blue color). The valence band 
states are QBs at the right boundary, and those of the 
conduction band are QBs at the left boundary. Though 
the band gap is very small, the low-frequency transitions 
between these two bands are drastically suppressed due 
to the very weak spatial overlap. This agrees with the 
corresponding optical spectra for £ = 0.144 V/A shown 
in Fig. |3] (a), which are almost the same as for £ = 0. 

Assume that the field is slowly increased to £ > £ m , 
so that the level anticrossing shown in Fig. Q] (f ) is adi- 
abatically passed, which requires that the passage time 
t p > h/AE m (cf. Landau-Zener theoryi^). This neces- 
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FIG. 3: Optical properties of nanofilms with N = 30 in elec- 
tric fields, (a) Below metallization threshold: real (blue) and 
imaginary (red) parts of permittivity e for fields £ < £ m . 
The dashed curves correspond to the zero field, while the 
solid curves are for £ — 0.144 V/A. (b) Above metalliza- 
tion threshold: Ree and Ime for £ = 0.152 V/A. Two low- 
frequency spectral peaks correspond to the intraband tran- 
sitions between the QBs. (c) The same as in the previous 
panel but in a stronger field £ = 0.238 V/A. (d) The same as 
(c) but for £ = 0.854 V/A. The three low-frequency absorp- 
tion peaks correspond to FKE, intraband transitions between 
WS states, and intraband transition between QBs, respec- 
tively, listed in the order of increasing transition frequency, 
(e) Combined oscillator strength of low-frequency transitions 
as a function of applied electric field £ . The green line is an 
analytical approximation - see text, (f) Real part of dielec- 
tric permittivity at zero frequency (blue) and imaginary part 
(red) of dielectric function averaged over low-frequency spec- 
trum (E < 3 eV) as functions of the applied field. Analytical 
approximations (green and orange lines) are described in text. 



sitates that the nanofilm is thin enough - cf. Eq. (|5|). 
Then the passage is adiabatic, the system persists in the 
ground state, and the electron population remains below 
the Fermi surface - see Fig. [2] (b). At each boundary 
of the film, there are QBs on the opposite sides of the 
Fermi surface, which significantly overlap in space. This 
allows for strong electron transitions at low frequencies 
as shown by the vertical double arrows. Note that the 
transition probability between two QBs rapidly decreases 
with the transition frequency ex a; -4 [see Eqs. (30) and 
(34) of Supplementary Material^ 1 -] . Correspondingly, the 
optical absorption is dramatically shifted to the red and 
infrared (ir) parts of the spectrum as displayed in Fig. 
[3] (b). In the red spectral region, Ree < 0, which is 
characteristic of metals. This is the metallization effect. 

With a further adiabatic increase of the field to £ = 
0.238 V/A, more levels cross the Fermi surface - see Fig. 
[2] (c). A stronger low- frequency absorption takes place 
[Fig. 0(c)] . This signifies a more developed metallization. 



For a very strong field case illustrated in Fig.[3](d), the 
Fermi surface separates levels originating from four ini- 
tial energy bands. The dominating transitions are those 
between WS states, and their frequency is ~ e£ a increas- 
ing linearly with £ . This leads to a general shift of the 
optical spectra to the blue - see Fig. [3] (d) where the 
absorption maximum is now at approximately 2 eV. 

To emphasize the metallization transition, we show in 
Fig. El (e) the oscillator strength / of the low frequency 
(below 3 eV) transitions as a function of £. Note that 
f/n sa 1 is characteristic of metals. This figure demon- 
strates a resemblance between the metallization and a 
quantum phase transition. Up to the metallization crit- 
ical field £ m — 0.148 V/A, the oscillator strength / is 
practically zero. After that, / increases in steps, each 
corresponding to a pair of QBs adiabatically crossing the 
Fermi surface. At the maximum, f/n ps 0.4 implying 
that ps 40% of the total number of electrons contributes 
to this metallic behavior. A further increase of £ leads to 
the WS states crossing the Fermi surface and the oscilla- 
tor strength decreasing [cf. Eq. (46) of the Supplemental 
Material 11 ] / oc £ ~ l as shown by the solid green line. 

Behavior of the permittivity across the metallization 
transition is displayed in Fig. [3] (f). Below £ m , e is rel- 
atively low. For £ > £ m , both Ree(0) and Ime increase 
in steps due to the QB states crossing the Fermi surface. 
After reaching the maximum, the main contribution to 
e shifts to the transitions between WS states originating 
from the same band (either valence or conduction) . This 
causes decrease in e as shown by the green and orange 
lines computed using Eqs. (47)- (48) of the Supplemen- 
tal Material 11 . As field increases, there are also transi- 
tions between the WS states originating from the differ- 
ent bands, leading to the sharp peaks in Fig. [3] (f). 

Concluding (see also Sec. VHI of Supplemental 
Material 11 ), we have predicted an effect of metallization 
in dielectric nanofilms induced by an adiabatically in- 
creasing applied field. The localized states crossing the 
Fermi surface cause optical absorption extending from 
very low (THz) frequencies over all optical region. In the 
near-ir and red spectral region, Ree < is predicted. 
This property is characteristic of metals and allow for 



a multitude of nanoplasmonic effects. The metalliza- 
tion cardinally differs from the Zener breakdown in bulk 
crystals, which is clear from much lower fields required 
(£ m ~ 0.1 V/A for L = 10 nm). In fact, the metalliza- 
tion is defined not by the field £ per se but by the total 
potential difference AC/ = £ m L — E g /e. Due to the re- 
quirement of adiabatic passage to the metallized state, 
the rise time t p of the applied electric field exponentially 
increases with L^J \xEgjh. For instance, for the consid- 
ered case (L = 10 nm, E g = 9 eV), t p > 10 ns, while for 
L = 3 nm the passage is much faster: t p > 1 fs. 

The manifestations of the metallization depend on the 
way the field is induced in the nanostructure. For the 
excitation by an optical or THz wave electric field, the 
metallization will cause high values of the permittivity 
and, consequently, bring about the plasmonic behavior 
of the system. This will lead, in particular, to screening 
of the external fields limiting the internal fields to £ ~ 
£ m . Note that £ m ~ 0.1 V/A corresponds to the wave 
intensity W ~ 10 11 — 10 12 W/cm 2 , which is well tolerated 
by nanostructured plasmonic metals - cf. Ref. Il7l In 
this case, the metallization effect is completely reversible. 
This will open up the field of nanoplasmonics to a variety 
of new dielectric and semiconductor nanosystems with a 
plethora of new phenomena possible. Among potential 
applications, is an ultrafast field-effect transistor where 
an ir or optical fs pulse controls a dielectric gate. 

In contrast, for dc- to microwave-frequency potential 
applied via electrodes, the external potential difference 
is fixed. Then the metallization will lead to and is a 
new mechanism of the dielectric breakdown, which is fun- 
damentally different from both the Zener and avalanche 
mechanisms. Such a situation is characteristic for the 
nanometric layers of the insulator in field-effect transis- 
tors and super-capacitors, with far-ranging technological 
ramifications for microelectronics and energy storage. 

We appreciate discussions with F. Krausz and R. Ern- 
storfer. This work was supported by grants from the 
Chemical Sciences, Biosciences and Geosciences Division 
of the BES Office of the US Department of Energy, a 
grant CHE-0507147 from NSF, and the US-Israel BSF. 



1 C. Zener, Proc. Royal Soc. A 145, 523 (1934). 

2 P. B. Corkum and F. Krausz, Nature Physics 3, 381 
(2007). 

3 F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163 (2009). 

4 T. Otobe, M. Yamagiwa, J. I. Iwata, K. Yabana, T. Nakat- 
sukasa, and G. F. Bertsch, Phys. Rev. B 77, 165104 (2008). 

5 G. H. Wannier, Phys. Rev. 117, 432 (1960). 

6 W. Shockley, Phys. Rev. Lett. 28, 349 (1972). 

7 F. Bloch, Z. Phys. A 52, 555 (1929). 

8 W. Franz, Z. Naturforschung A 13, 484 (1958). 

9 L. Keldysh, J. Experimentl. Theor. Phys. 34, 5, 1138-1141 
(1958); Translation: Sov. Phys. JETP 7, 788 (1958). 

10 D. A. B. Miller, D. S. Chemla, and S. Schmitt-Rink, Phys. 



Rev. B 33, 6976 (1986). 

11 See supplementary material at arxiv.org/abs/1007.2366 
(2010). 

12 D. M. Goodmanson, Am. J. Phys. 68, 866 (2000). 

13 E. O. Kane, J. Phys. Chem. Solids 12, 181 (1959). 

14 K. B. McAfee, E. J. Ryder, W. Shockley, and M. Sparks, 
Phys. Rev. 83, 650 (1951). 

15 S. Glutsch, Phys. Rev. B 69, 235317 (2004). 

16 C. Zener, Proc. Royal Soc. A 137, 696 (1932). 

17 S. Kim, J. H. Jin, Y. J. Kim, I. Y. Park, Y. Kim, and S. W. 
Kim, Nature 453, 757 (2008). 



